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1.  Introduction 

Measurement  of  vehicle  dynamics  during  flight  on  test  ranges 
has  been  approached  using  both  ground-based  and  onboard  instrumentation. 
Onboard  platforms  and  accelerometers  require  telemetry  systems  or 
recorders  resulting  in  relatively  expensive  instrumentation  which  must 
be  considered  expendable  in  most  tests.  Onboard  solar  aspect  sensors 
have  been  developed  for  spinning  missiles  which  can  be  used  to  infer 
yawing,  pitching,  and  rolling  motion  of  flight  vehicles.  The  most 
successful  of  these  is  the  yawsonde  developed  at  the  Ballistic  Research 
1 2 

Laboratories  * . Many  flight  systems  are  tested  using  only  ground-based 
instrumentation  employing  primarily  photographic  techniques.  Data 
reduction  in  these  cases  is  essentially  all  manual  with  attendant  high 
cost  and  relatively  long  data  reduction  time. 

Two  approaches  for  the  determination  of  position  and  attitude  of 
flight  vehicles  are  described  in  this  report.  These  systems  involve 
ground-based  laser  transmitter/detector  stations  and  expendable  retro- 
reflecting  elements  located  on  the  vehicle.  When  compared  to  current 
ground-based  photographic  techniques  for  determining  both  position  and 
attitude,  the  proposed  systems  offer  several  advantages  including: 
lower  recurring  costs;  improved  accuracy;  automated  data  reduction 
and  applicability  in  low  ambient  light  situations,  and  in  situations 
where  the  trajectory  of  the  vehicle  is  not  well  established  prior  to 
flight. 

This  report  describes  the  two  concepts  and  presents  some  limited 
results  that  were  obtained  from  system  analysis  studies. 


2.  The  Three  Station  System 
a . System  Concept 

This  system  includes  three  ground  based  transmitter/detec- 
tor  tracking  stations  - each  incorporating  one  pulsed  and  one  continuous 
wave  (CW)  laser  of  different  frequencies.  Onboard  the  vehicle,  two 
different  types  of  retroref lecting  arrays  are  required. 


Mermagen,  W.  H.  and  Clary,  W.  H. , The  Design  of  a Second  Generation 
Yawsonde,  Ballistic  Research  Laboratories,  Aberdeen  Proving  Ground, 
Maryland,  BRL  Memorandum  Report  No.  2368,  1974  (Unclassified). 
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One  type  of  array  is  composed  of  eithei  conventional  corner  cubes, 
reflective  tapes  and/or  paints  which  have  the  property  of  retroref lect- 
ing  a portion  of  the  collimated  incident  beam  back  parallel  to  itself 
regardless  of  the  orientation  of  the  reflective  surface.  This  array 
forms  a retroref lecting  band  located  on  the  perimeter  of  the  vehicle 
body  at  one  axial  position.  Illumination  of  and  reflection  by  this 
retroref lecting  band  will  then  form  a conventional  laser  radar  tracking 

„ 3,4 

system 

A series  of  roof  type  prisms  as  shown  in  Figure  1 form  a second 
array  which  is  also  mounted  on  the  vehicle's  surface.  This  particular 
reflector  array  will  be  designated  here  as  the  single  plane  corner 
reflector.  A plane,  which  passes  through  the  center  of  all  the  reflect 
ing  surfaces  of  the  roof  type  array  and  is  also  normal  to  all  these 
surfaces,  will  be  called  the  retroref lection  plane.  Collimated  light 
incident  on  the  single  plane  corner  reflector  and  contained  in  its 
retroref lection  plane  (£  = 90°)  is  reflected  back  to  the  source.  The 
single  plane  corner  reflector  is  mounted  on  the  surface  of  the  vehicle 
in  such  a manner  that  its  retroref lection  plane  contains  the  roll  axis 
of  the  vehicle.  An  illustration  of  an  instrumented  vehicle  is  pro- 
vided in  Figure  2. 

As  the  vehicle  flies  downrange,  it  is  tracked  with  the  laser  radar 
which  positions  a CW  laser  to  provide  a continuous  illumination  of  the 
vehicle  and  also  provides  the  information  for  determining  the  vehicle 
position  as  a function  of  time.  During  each  revolution  of  a spinning 
vehicle,  a CW  laser  pulse  is  returned  to  every  tracking  station.  The 
time  interval  between  the  pulses  returned  to  three  separate  tracking 
stations  provides  enough  data  for  determining  the  vehicle  attitude. 

A mathematical  description  of  this  system  is  presented  in  the  next  sec- 
tion. 


It  should  be  noted  here  that  the  missile  position  can  be  determined 
using  triangulation  from  two  or  more  tracking  stations.  Range  determi- 
nation with  the  use  of  a pulse  laser  is  therefore  not  required  for  a 
multistation  system.  A single  CW  laser  could  be  used  for  both  attitude 
and  position  determination  if  discrimination  between  the  two  types  of 
reflections  were  possible. 
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b.  Mathematical  Description  of  the  System 

This  analysis  assumes  constant  values  for  the  vehicle 
position,  pitch,  yaw,  and  roll  rate  during  the  time  it  takes  for  a CW 
pulse  to  be  received  at  each  of  the  three  tracking  stations.  The  impact 
of  these  assumptions  will  be  treated  as  an  error  in  the  subsequent  error 
analysis.  A vehicle  which  is  fixed  in  space,  rotating  at  a frequency 
w,  equipped  with  a previously  described  single  plane  corner  reflector 
and  illuminated  with  a CW  laser  will  retroreflect  the  incident  CW  beam 
to  a detector  co-located  with  the  CW  laser  source  and  thereby  generate 
a periodic  signal  of  frequency  w. 
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SINGLE  PLANE 
CORNER  REFLECTOR 


Figure  2.  Illustration  ot  ''ehicle  showing  retroref lectors . 


The  mat!  ematical  description  may  be  facilitated  by  defining  three 
coordinate  systems.  An  earth  fixed  Cartesian  system  (Figure  3)  is 
defined  with  the  following  characteristics: 

1)  The  origin  is  located  at  the  launch  site. 

2)  y positive  is  in  the  downrange  direction. 

3)  x is  crossrange. 

4)  z is  positive  in  the  vertical  upward  direction. 

5)  (x,y,z)  forms  a right  hand  system. 

A vehicle  based  coordinate  system  (Figures  3 and  4)  is  defined 
such  that: 

1)  The  origin  is  located  at  the  vehicle  center  of  gravity. 

2)  coincides  with  vehicle  roll  axis  with  the  positive  direction 
toward  the  nose. 

3)  r,  is  perpendicular  to  u;  and  parallel  to  the  x-y  plane  of  the 
earth-fixed  system. 

4)  I is  perpendicular  to  and  w such  that  (’,,u>,'l  is  a right- 
hand  system. 

An  intermediate  system  (x1 , y' , z')  which  is  just  a direct  translation 
of  the  earth  fixed  system  to  the  vehicle  center  of  gravity  is  also 
defined  as  shown  in  Figure  4. 

One  approach  to  analysis  of  the  system  involves  a description  of 
the  locations  of  the  ground  transmitter/detector  stations  in  terms  of 
the  vehicle  based  (-,,10,1)  coordinates.  Let  R.  be  the  position  vector 


of  the  i- 


.th 


ground  tracking  station  in  the  earth  fixed  system,  R be 

m 


the  position  vector  of  the  vehicle  in  the  earth  fixed  system,  and  R.' 

th  ^ 

be  the  position  vector  of  the  i—  ground  tracking  system  in  the  x1 , y' 

and  z'  system.  That  is, 


R. ' = R.  - R 
i i m 

.th 


(1) 


The  position  vector  of  the  i — ground  tracking  station  in  the 

vehicle  based  coordinate  system  (n.,  w.,E  ) will  be  denoted  as  . The. 

i 1 l i 

relationship  between  y.  and  R.1  is  given  by 


= (a)  R. 

i l 


(2) 
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, y1  V coordinate  system. 

Figure  4.  X , > 

relative  to  earuti  ^ 

is  *l»  — «*  ^ «a„s£o™ation  t» 

® ^ e,oatiros  . <*> 

(b  ) (x  - x ) - sinC&2)  " yrn 
= cos(&2)  i*i  m , _ Y ) (5) 

lb)l%  - X ) + cos(82>  C0S(61)  &i  “ 

, = sin(%)  cos&J  C i m 

i 4 

+ sin(-^)  ^zi  Zm^ 
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| = -sin(5  ) sin(6  ) (x  - x ) - cos(£>  ) sin(6  ) (y.  - y ) 

1 l 1 1 tr;  4 i 1 ID 

(6) 

+ cos(&, ) (z.  - z ) 
i i m 

which  gives  a relation  between  the  vehicle's  measured  position  (x.  - x , 

l m 

y.  - y^,  z^  - z^)  and  the  5 unknown  parameters:  pitch  (6^,  yaw  (&2), 

and  the  station  coordinates  w.  , r,^). 


r 


The  two  additional  relationships  needed  to  close  this  system  of 
equations  can  be  developed  by  determining  the  angle  that  the  retrore- 
f lection  plane  makes  with  the  q coordinate  which  is  denoted  as  0 in 
) .gure  5.  Since  the  roll  axis  w is  contained  in  the  retroref lection 
plane,  this  angle  can  be  expressed  in  terms  of  the  vehicle  based  coordi- 
nates q and  i.  Further,  the  angle  swept  out  by  the  retroref  lection 
plane  as  it  moves  between  two  ground  stations (i  and  j)  can  be  expressed 
in  terms  of  these  same  two  coordinates  as  indicated  in  Figure  5. 


r.‘ 


|i  . i 


m 


<t>.  - ❖ . = tan  | ) - tan  M (7) 

1 J \\l  Y»j  / 

This  can  be  equated  to  the  product  of  the  roll  rate  (w)  and  the  measured 
time  lapse  (At  ) between  the  pulses  received  at  two  of  the  ground 

stations  when  the  retroref lection  plane  passes  through  them. 

<t> . - $ . = uAt . . (8) 

i J JO 

Any  two  of  the  three  possible  time  differences  which  can  be  formed  in 
this  manner  (At » At^  and  At^)  cou^  be  used  to  close  the  above  set 

of  equations  (4,  5,  and  6).  Choosing  At^  anc*  At^  and  making  the 

appropriate  substitutions  reduces  this  set  of  equations  to  the  following 
two  equations: 

1 -slnC,)  Sin(f  ) (x,  - x > - cos{t  ) sln(0  (y,  - y ) + cos(i-,)  (z,  - z ) 

.'t.,  - i arccan  A—HL 2 J 2 E- ! 1 E_ 

12  ^ cos(6,)  (x,  - X ) - sin(b  ) (y,  - y ) 

| * t m 4 £■  tn 


-sin(f>  ) sin(^  ) (x  - x ) - cos(i  ) sln(£.)  (y.  - y ) + cos(r 
4 i in  4 i 1 rn 


cos C.)  (x  - x ) - sln(b.)  (y  - y ) 
4 i rn  4 l m 


) (z.  - z )") 
1 m 


-sln(t  ) sln(1-  ) (x.  - x ) - cos(b„)  sln(b  ) (y  - y ) + cos(t.  ) (z  - z ) 
410  m 4 l J tn  i Otn 

cos(b  ) (x,  - x ) - stn(i,)  y.  - y ) 

4 0 m 4 J m 

ln(b  ) sinC  ) (x  - x > - cos(r  ) stn('  ,)  (y.  - y ) + cosC  ) (z  - z ) 

4 llm  4 l l m 1 i rn  \ 

cos(V  (Xj  - xra)  - sin(  2)  (yt  - yn>  J j 


BSMfewS B 


* 'i'r <'jMft»  v 


!B 


M 


i 


a 


Oi*?» 


■ Vts 
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w can  be  determined  from  the  time  interval  between  successive 
pulses  received  at  a ground  station  and  the  same  pulse  information  is 


used  to  calculate  the  time  differences  At^  and  At^. 


The  vehicle's 


position  relative  to  the  ground  stations  (x.  -x,y.  -y,z.  -z) 
r 1 m Jx  m 1 m 

is  measured  by  the  tracking  part  of  the  system.  The  pitch  (6^)  and 
yaw  (6^)  °f  the  vehicle  can  therefore  be  determined  by  the  simultaneous 
solution  of  the  previously  noted  equations. 

These  two  equations  are  transcendental  and  may  possibly  have  mul- 
tisolutions but  this  does  not  seem  to  be  the  case  in  the  investigations 
carried  out  to  date. 

c.  Ground  Station  Position  Optimization 

Three  major  considerations  in  the  optimization  of  ground 
station  locations  for  the  laser  systems  were  used: 

1)  In  view  of  the  high  accelerations  and  resultant  large  veloci- 
ties associated  with  rockets  and  missiles,  the  stations  must  be  located 
so  that  the  required  angular  tracking  rate  of  the  azimuth/elevation 
mount  will  be  within  the  mount's  capabilities. 

2)  The  efficiency  of  a single  plane  corner  reflector  (roof  prism) 
decreases  rapidly  as  the  angle  of  incidence  (0  in  Figure  1)  diverges 
from  90°.  Thus,  the  intensity  of  the  reflected  beam  decreases  from 
approximately  1007,,  of  the  incident  beam  at  normal  incidence  (0  = 90°) 
to  zero  at  0 = 45°.  This  dictates  ground  station  locations  which  pro- 
vide working  angles  relative  to  the  missile  axis  (if  the  reflecting 
surfaces  are  arranged  symmetrically  with  respect  to  the  missile  axis) 
ranging  from  60°  to  120°  in  order  to  have  the  reflected  beam  intensity 
greater  than  approximately  507»  of  the  incident  beam. 

3)  The  third  constraint  arises  from  the  fact  that  the  equation 

for  At . . in  terms  of  pitch  and  yaw  cannot  be  solved  explicitly  for 
" J 

pitch  and  yaw.  Simultaneous  numerical  solution  of  the  two  equations 
is  facilitated  if  the  ground  stations  are  arranged  so  that  each  of 
the  partial  derivatives  dAt„/d5^  and  ^At^VdS^  is  either  positive 

or  negative  over  the  entire  dynamic  range  of  the  vehicle  flight.  These 
pitch  and  yaw  sensitivity  expressions  are  developed  in  the  next  session. 

To  check  the  static  model  developed  above  in  a dynamic  situation, 
a six  degree  of  freedom  computer  simulation  of  a particular  rocket 
(ARROW)  was  modified  to  include  subroutines  that  provide  return  pulse 
data  and  sensitivity  values  for  a simulated  laser  attitude  system. 

After  a number  of  simulations  it  was  determined  that  a static  missile 
position  could  be  used  to  determine  the  sign  of  the  sensitivities 
developed  from  the  At„  equations.  This  simplification  results  in  a 
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considerable  reduction  in  computer  time  since  the  trajectory  simulation 
does  not  have  to  be  exercised  for  every  ground  station  geometry,  and 
further,  it  allows  a more  systematic  study  of  the  interaction  of  missile 
position  and  ground  station  geometry  with  respect  to  the  various  sensi- 
tivities . 

A computer  program  was  developed  which  produced  the  various  sensi- 
tivities for  fixed  missile  position  and  attitude  and  a particular  ground 
station  geometry.  The  missile  position  was  varied  over  the  grid: 

x =0. 
m 

ym  = 10,  100,  200,  300,  400,  500,  1000,  1500,  2000  (ft). 

zm  = 10  (ft)- 

m m 


The  attitude  was  varied  over  the  grid: 
6^  = 0 to  10  in  increments  of  1°. 
5„  = -5  to  5°  in  increments  of  1° . 


Roll  rate  was  held  constant  at  3020°/sec.  Ground  station  positions 
were  varied  over  the  domain: 


x.  = -3000  to  3000  ft  in  increments  of  1000  ft. 

i 

y.  = 0 to  3000  ft  in  increments  of  1000  ft. 

i 

z . = 0. 
i 


Pitch  and  yaw  sensitivities  were  calculated  for  121  pairs  of  pitch  and 
yaw  values  for  each  combination  of  ground  station  locations. 


For  the  yaw  equation  (i.e.,  the  equation  in  which  pitch  is  held 
constant  in  an  iterative  solution),  the  program  printed  out  the  arith- 
metic mean,  the  maximum,  and  the  minimum  values  for  pitch  and  yaw  sensi- 
tivities if  the  yaw  sensitivities  were  all  either  greater  than  or  less 
than  zero.  If  values  of  both  signs  occurred,  the  next  combination  of 
tracking  stations  was  checked  without  providing  a computer  printout  of 
the  case  which  failed  criteria  3.  Printout  was  also  suppressed  if  the 
missile  pointed  directly  at  a ground  station  due  to  a singularity  in 
the  sensitivity  equations  for  this  case.  For  the  pitch  equation,  a 
similar  output  was  generated. 
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This  study  indicated  that  a wide  range  of  ground  station  positions 
were  available  for  attitude  determination  but  no  single  combination  of 
three  ground  station  positions  was  found  which  would  provide  coverage 
over  the  complete  dynamic  range  of  the  missile  as  specified.  In  retro- 
spect, a negative  pitch  attitude  during  the  first  2000  feet  of  flight 
appeared  unrealistic,  and  since  it  was  determined  that  these  values  were 
excluding  a large  number  of  possible  station  locations,  the  investigation 
was  repeated  constraining  5^  range  from  0 to  10°  in  increments  of  1°. 

The  results  from  the  previous  investigation  indicated  that  the  location 
(3000,  0,  0)  would  be  a reasonable  position  for  the  first  ground  station. 
Using  this  location  for  station  one,  86  combinations  for  the  other 
two  ground  station  positions  providing  positive  results  in  the  sense 
of  criteria  3 were  developed.  An  examination  of  the  mean  of  the  calcu- 
lated yaw  sensitivities  indicated  that  maximum  values  occurred  when 
one  ground  station  was  located  along  the  line  y = 3000  feet.  Figure  6 
presents  a plot  of  yaw  sensitivity  as  a function  of  a ground  station 
position  along  the  line  y = 3000  where  the  other  two  stations  are  located 
at  (3000,  0,  0)  and  (3000,  2000,  0).  In  addition,  plots  of  the  pitch 
and  yaw  sensitivity  envelopes  as  a function  of  the  downrange  missile 
position  for  two  other  combinations  of  tracking  station  positions  are 
shown  in  Figures  7 and  8. 

Considering  the  second  constraint  in  the  light  of  the  preceding 
results,  it  is  apparent  that  at  least  two  retroref lecting  prisms  with 
different  orientations  of  the  single  plane  corner  reflectors  will  be 
required  on  a vehicle  in  order  to  provide  coverage  downrange  for  2000 
feet.  Although  actual  constraints  on  angular  tracking  rate  (criteria  1) 
were  not  investigated,  the  angular  rate  for  a third  station  located 
anywhere  along  the  line  y = 3000  would  be  less  than  that  required  of 
stations  1 and  2 which  apparently  can  be  handled  with  off-the-shelf 
equipment. 

d.  Error  Analysis 


In  a dynamic  situation  where  a retroref lector  equipped 
rocket  flies  a trajectory  such  that  the  laser  pulses  are  returned  to 
a number  of  ground  stations , the  time  interval  between  reception  of 
pulses  at  any  two  ground  stations  is  influenced  by  changes  in  the  missile 
position  and  attitude  between  the  times  of  reception  of  the  pulses  by 
the  two  stations,  as  well  as  changes  in  the  roll  rate  during  the  time 
interval. 


In  the  analysis  that  was  presented  in  Section  2.b.  for  determining 
the  missile  pitch  and  yaw  from  the  points  in  time  when  pulses  were 
reflected  bade  to  one  of  the  three  gound  stations,  it  was  assumed  the 
missile  position,  attitude  and  roll  rate  did  not  change  during  the 
time  intervals  associated  with  the  reception  of  pulses  at  the  three 
ground  stations.  In  the  simplified  analysis  to  follow,  the  observed 
time  interval  At.,  between  any  two  station  pulse  receptions  is  treated 
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Figure  6.  y/J^  possible  third  station  locations  based  on  criteria  3. 

as  a nominal  time  interval  At„  plus  a number  of  perturbation  terms. 

At.,  is  defined  as  the  time  which  would  be  observed  if  the  missile  were 
ij 

fixed  in  space  at  a position  midway  between  the  positions  associated  with 
the  missile  at  the  times  of  reception  for  the  two  pulses;  possessed  a 
pitch  and  yaw  equal  to  the  mean  of  the  values  associated  with  the 
vehicle  at  the  times  of  reception  for  the  two  pulses;  and  had  a roll 
rate  equal  to  the  mean  of  the  values  associated  with  the  vehicle  at 
the  times  of  reception  for  the  two  pulses.  The  perturbation  terms  are 
developed  assuming  that  the  influencing  parameters  (variables)  are 
independent  and  further,  that  the  influence  can  be  represented  as  a 
linear  variation  arourJ  the  nominal  value;  in  other  words,  as  the  first 
term  in  the  Taylor  Series  expansion.  Mathematically  this  may  be  repre- 
sented as 


15 


At . . = At . . + 
1J  1J 


i Rh-v 

K-l  L.  t=t 


UJl  . , 

where  -g^— *•  is  the  sensitivity  of  the  time  interval  At. . to  variations 

k ^ 

of  parameters  <t>,  which  includes:  4>,  = 5,  , <t>n  = 5. , <t>_  = to,  <t>  = x , 

k 1 1 l i 3 4 m 

= y > = 2 > anc^  = I.  <J>,  is  the  nominal  value  of  the  parameter 

o m o m / k r 

associated  with  At„.  1 is  an  instrumentation  parameter  which  incorpo- 
rates all  of  the  system  uncertainties  such  as  clock  accuracy. 

The  maximum  error  that  vehicle  dynamics  and  instrumentation  intro- 
duce in  the  determination  of  5^  will  be  investigated  in  the  following 


dAt . . " „ 9 r-v  | At . . 
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The  various  sensitivities  are  developed  in  the  following  subsections. 

(1)  Pitch  Sensitivity.  The  pitch  sensitivity 

o5j 

can  be  obtained  by  using  Equations  (7)  and  (8). 


- arctan  — 
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the  pitch  sensitivity  is 
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Since  r,^  is  independent  of  6^  [Equation  (4)]  one  obtains 
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(2)  Yaw  Sensitivity, 


^ (At) 

d5„ 


. Approaching  the 


nation  of  yaw  sensitivity  in  the  same  manner,  one  obtains 


m 


l 

w 


= A 


. dq.  dq^ 

d&2  db2  ^i  db2  ^i  d&2 


r,2  + t2 

'j  J 


,2+  ?2 

'l  1 


where 


dt . 

T£T  = -cos (6,)  sin (5  ) (x.  - x ) + sin(6  ) sin(&  ) (y.  - 
oo2  L l i m 2 1 i 


dS„ 


l _ 


- -sin(6  ) (x  - x )-cos(5  ) (y.  - y ) 
z 1 m ~ 


Z i m 


19 




"Wt»:  ~ 


(16) 


) 

(17) 


determi- 


(18) 


(19) 

y ) (20) 

m 

(21) 


'iff 

1% 


tavs meZg&HA. 


S-y  ^ vtf 


(3)  Altitude  Sensitivity,  ^LQl.  t The  altitude  sensi- 
tl  /ity  is  given  by  OZm 
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and  since  *•,.  is  not  a function  of  z 
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tivity  is  simply 


(4)  Roll  Rate  Sensitivity,  The  roll  rate  sensi- 


d(At..)  . 

n _ -1 


arctan 


/ii\  „ fh\] 

\'j/  WJ 


^At l 


given  by 


(5)  Range  Sensitivity,  $ y • The  range  sensitivity  is 
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(6)  Crossrange  Sensitivity. 
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In  order  to  proceed  further  in  a practical  sense,  some  estimate 
of  the  values  of  (0^  - 0^)  must  be  determined  to  provide  numerical 

values  for  the  terms  in  the  summation,  as  well  as  At,  .,  and  At, . in 

ij 

Equation  (12).  The  approach  taken  in  this  investigation  has  been  to 
choose  a representative  vehicle  (ARROW),  and  exercise  a six-degree-of- 
freedom  computer  simulation  to  provide  some  representative  values. 

Assuming  three  ground  stations  with  coordinates  (3000,  0,  0), 
(3000,  2000,  0),  and  (1000,  3000,  0),  the  results  of  the  simulation 
of  an  ARROW  flight  indicate  that  the  - 4^  are  of  the  orders  of 

magnitude  indicated  in  Table  1.  The  various  sensitivities  are  also 
ca:;ulated  in  the  simulation  program.  Using  the  calculated  sensitivi- 
ties and  the  4>,  - 4>  determined  from  the  simulation,  an  estimate  of 

K K 

the  error  arising  from  vehicle  dynamics  may  be  obtained.  In  the 
example  used  (ARROW),  the  dynamics  of  the  vehicle  introduced  an  error 

-4 

in  pitch  of  approximately  8 X 10  degrees  and  in  yaw  of  approximately 
_2 

6 X 10  degrees. 
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TABLE  1.  ERROR  ANALYSIS  SUMMARY 


♦k 

Pitch  System 

Yaw  System 

|jatij/o»kl 

l*k  - 

l^ij/^kl 

K - \\ 

Roll  rate  (°/ sec) 

6.702E-8 

1.118 

3.796E-8 

4.160 

X 

m 

9.625E-9 

7.355E-5 

6.639E-9 

2.737E-4 

y 

m 

, ft 

2.201E-9 

5.526E-2 

1.462E-9 

2.057E-1 

z 

8.364E-8 

1.176E-3 

2.370E-4 

8.752E-3 

m 

Pi tch  \ 

3.457E-4 

1.730E-5 

Yaw 

j deg 

8.804E-7 



b.  Conclusions 


This  study  indicates  that  a three  ground  station  system 
employing  retroref lecting  prisms  (single  plane  corner  reflectors)  can 
provide  a useful  range  instrumentation  technique  for  vehicle  attitude 
determination.  Yaw  sensitivity  and  therefore  yaw  accuracy  will  be 
.ess  thar.  pitch  accuracy.  Furthermore,  for  the  particular  combination 
of  station  locations  investigated,  the  system  sensitivity  is  a func- 
tion of  downrange  position  with  an  indication  of  severe  degradation  in 
> aw  sensitivity  for  some  combinations  of  flight  parameters  near  launch. 
It  is  possible  that  a portable  third  station  could  be  located  at  an 
optimum  position  for  maximizing  yaw  sensitivity  for  each  flight  test. 
However,  the  three  station  system  investigated  combined  with  optical 
lever  instrumentation  for  the  near  launch  phqse  of  flight  may  provide 
sufficiently  accurate  flight  test  data  through  the  first  2000  feet  of 
flight. 

A number  of  practical  problems  associated  with  the  overall  system 
accuracy  have  not  been  addressed.  Specifically,  tracking  accuracy,  CW 
laser  transmitter  and  detector  characteristics,  tolerance  of  optical 
components,  fabrication  tolerances,  preflight  calibration,  deformation 
of  the  reflectors  due  to  acceleration,  and  thermal  stress  were  not 
considered.  None  of  the  above  considerations  appear  to  present  insur- 
mountaole  problems;  however,  they  all  will  contribute  to  degradation  of 
the  system  performance  and  therefore,  warrant  further  consideration. 


3.  The  Two  Station  System 

The  requirement  of  a third  station  can  be  eliminated  using  a 
two  station  concept  which  requires  two  single  plane  corner  reflectors 
on  the  vehicle  mounted  in  such  a fashion  that  the  planes  of  tie  corner 
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reflectors  (defined  previously)  are  skewed  relative  to  one  another.  As 
the  vehicle  is  tracked  downrange  two  pulses  are  returned  to  each  track- 
ing station  for  each  revolution  of  the  vehicle.  The  time  intervals 
between  the  various  pulses  again  provide  the  data  to  determine  pitch, 
yaw,  and  roll  attitude. 

a.  System  Concept 

The  geometrical  arrangement  of  the  optical  components  is 
considered  as  shown  in  Figure  9.  The  two  single  plane  corner  reflectors 
are  located  at  the  same  axial  position  on  the  vehicle  but  are  separated 
by  a polar  angle  equal  to  2p.  Further,  they  are  inclined  relative  to 
the  roll  axis  by  angles  7 ^ and  y0 . As  a result  the  planes  of  the  two 

corner  reflectors  are  skewed  with  respect  to  both  the  roll  axis  and  each 
other  as  shown  in  Figure  10.  When  the  vehicle  rotates  to  a position 
so  that  a retroref lection  plane  passes  through  a transmitter/detector 
ground  station,  a pulse  is  retroref lected  to  the  ground  station.  Con- 
sider a vehicle  fixed  in  space  and  rotating  at  a frequency  w.  Referring 
to  Figure  9,  it  may  be  seen  that  the  time  interval  between  retroref lected 
pulses  are  given  by 


2jt(iAt^2  = 2 P + 0^  + 


(33) 


Note  from  Figure  10  that 


tan  a = - , 


where  a is  called  the  laser  aspect  angle  and 

C. 

tan  ri  ' ~t 


so  that 


(C  ) (tan  7.) 

0.  = arcsin  = arcsin  —rr C" 

i (R)  (tan  a) 


Equation  (33)  may  be  rewritten  as 


2nw^t^2  - 2(3  = arcsin 


tan  7^ 


tan  a 


+ arcsin 


tan  yr 


tan  cr 


(34) 
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Figure  9.  Corner  reflector  geometry  relative  to  the  vehicle. 
Taking  the  sin  of  both  sides,  one  obtains 

''tan  yr 


sin  {2  Grw^t  -p)}  = — cos  | 

12  tan  a I 


arcsm 


tan  0 


tan  7 £ 


tan  cr 


co 


15  {arcsin  (tst)}  • 


(35) 


Noting  **he  geometry  show,  in  Figure  11,  the  following  expression  is 
obtained  : 


arcs 


in 


/tan  7 A | (tan2  a - tan2  7i)1^2 


cos 


(tan^o  - tan^7j)^  ^ 

Figure  11.  Geometry  for  substituion. 


so  that  (35)  may  be  expressed  as 

tan  7 (tan2  a - tan2  7.) 1/2 

sin  {2(-r^t  - (3)}  = ^ 2_ 

12  tan  a tan  a 


2 2 1/2 
tan  7^  (tan  a - tan  7^) 


tan  a 


tan  a 


n2cr  sin  {2(iu£t12  - P)}  = tan  ?1  (tan 2o  - tan2  72)1/2 

2 2 1/2 
+ tan  72  (tan  a - tan  7^)  ' 


Now  taking  the  cos  of  both  sides  of  (34)  one  has 


cos  {2(nu&t:i2  - p)}  = 


(tan^  a - tan2  7^ 


i2  (tan2  a - tan2  72)1//2  - tan  7? 

~2  “ 

tan  a 


This  may  be  solved  for  (tan2  a - tan2  7 )^2  to  yield 


2 2 1/2 
(tan  a - tan  7^)  = 


(2(ru&t^2  - P)}  tan2cr  + tan  7 1 

2 ~ 2 .1/2 

(tan  a - tan  70) 


cos 


When  substituted  into  Equation  (38)  one  obtains 

sin  - P)}  tan2a  = tan  7 (tan2a  - tan2  7^) 


1/2 


[cos 


+ tan  y. 


{2  (nu&t^  - P)1  tan^a  + tan  7^  tan  72 


2 2 1/2 

(tan  a - tan  72) 


2 2 1/2 

Multiply  both  sides  by  (tan  a - tan  72)  gives 

2 2 1 /2  \ 2 2 

(tan  a - tan  y^)  ' sin(2(nojAt  - P)}  tan  a = tan  7 tan  a 

+ tan  7 2 cos(2(noi^t^2  ~ P))  tan2<j 


Dividing  through  by  sin  ^ lUJLiL,  « 
yield 


(41) 


(42) 


tan  a - tan  7 2 = 


2 

(2(-tui^t^2  “ P))  tan  a and  squaring  both  sides 
{tan  y1  + tan  7 2 cosUtjiu&t^  - P))}2 


sin2{2(nw^t12  -P)) 


. (43) 
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This  may  be  solved  for  tan  0: 

2 I 2 

tan  a = {tan  7^  + 2 tan  7^  cos(2  (rtwAt^  - P))  tan  7 2 


+ tan  72  cos 


2{2(nw^t12  - P)1  + tan2  7 2 sin2  {2(«u&t12  " P)l] 

(44) 


sin2  (2(«u^t12  - P)) 


Alternately 


2 I 2 

tan  cr  = { tan1  7^  + 2 tan  7^  tan  7 ^ cos{2(nuAt^2 


P)}  + tan2  7 2J 


(45) 


sin 


2{2(nw^t12  - P)} 
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Take  the  square  root  of 
yielding  tan  a = tan  7^ 


tan  7 

both  sides,  let  = A,  and  2nw££  2f3  = (X 

’ tan  7^ 

2 1 /? 

(1  + A + 2A  cos  0:)  '/sin  Q.  (46) 


It  is  apparent  that  Equation  (46)  provides  a relationship  for  the 
laser  aspect  angle  (a)  which  is  the  angle  between  the  roll  axis  of  the 
vehicle  and  the  vector  joining  the  tracking  station  and  the  vehicle. 
This  relationship  is  in  terms  of  fixed  geometric  parameters 

and  two  parameters  which  are  measured,  i.e.,  u (determined  from  the 
time  interval  between  pulses  returned  to  one  ground  station  from  one 
corner  reflector)  and  Lt  (determined  from  the  time  interval  between 
pulses  from  the  two  different  corner  reflectors). 

b.  Data  Reduction 


For  simplicity,  a configuration  will  be  considered  which 
has  one  reflector  parallel  to  the  roll  axis,  a second  reflector  skewed 
an  angle  y to  the  roll  axis,  and  two  laser  tracking  stations.  Two 
types  of  return  pulses  are  received  at  the  stations.  The  time  interval 
for  return  pulses  from  the  unskewed  single  plane  corner  reflectors  to 
the  two  ground  stations  is  given  by  Equation  (9)  which  is  rewritten 
as 


ut 
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-sln(&2)  sln(6j)  (x2  - x^)  - cos(t>2)  sln(&t>  y2  - yj  + cos(6j)  (z2  - z^) 


cos  ( . 2 ) u2 


X ) - sln(&,)  (y  - y ) 
m 2 2 m 


-sin(6,)  sin(6  ) (x  - x ) - cos(&  ) sln(&.)  (y,  - y ) + cos(&.)  (z,  - z ) 
2 1 l m 2 1 l m 1 1 in 


cos (b9)  (x  - x ) 
2 l tn 


sln(&  ) (y  - y ) 
2 1m 


(47) 


Referring  to  Figure  5 it  may  be  seen  that  the  aspect  angles  can  be 
expressed  in  terms  of  qto|  coordinates  such  as 


tan  cr^ 


w 


1 


(48) 


Substituting  in  Equation  (15)  and  using  the  expressions  for 
given  previously  (Equations  (4),  (5),  and  (6)),  one  obtains 


cosb  (x.  - x ) + cosb.  cosb,(y,  - y ) ♦ sinb,(?,  - t ) 

t 1 . n, £ i 1 . -3. 1 _ 1,,  ,a_ 


£ — ft £ — 2 3 T77T  • ('AQ  ) 

|(co*&2(xl  - xo)  - Sin  2(y1  - ya>i  + Usinv-2  slnb^Xj  - xq)  - cosb2  tinb^yj  - yo>  + cost^z  - tji  } r ^ ' 
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It  is  apparent  that  Equations  (47)  and  (49)  form  a pair  of  simul- 
taneous equations  in  two  unknowns , i.e. , 5^  and  As  in  the  three 

station  system  a scheme  such  as  the  Newton-Raphson  technique  is  required 
to  determine  5^  and  82*  Recent  studies  indicate  that  multivalued  solu- 
tions exist  for  this  pair  of  equations.  Although  it  appears  that  physical 
constraints  may  be  used  to  eliminate  at  least  one  of  the  possibilities, 
a definitive  statement  in  this  regard  cannot  be  made  at  this  time.  It 
should  be  noted  that  it  may  be  possible  to  follow  a unique  solution 
curve  if  continuous  data  from  launch  is  available.  This  area  is  currently 
under  investigation. 

c.  Sensitivity  and  Accuracy 

A detailed  investigation  of  sensitivity  and  accuracy 
particularly  with  reference  to  vehicle  dynamics  has  not  been  completed 
for  the  two  station  concept.  A study  of  the  errors  associated  with  a 
determination  of  a using  Equation  (46)  is  presented.  Rewriting 
Equation  (33),  one  obtains 


At12  = (2P  + Bj  + B2}  ^ 


(50) 


or  alternately 
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f 

/tan  7 \ 

/tan  7 \ | 
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1 
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(51) 


The  sensitivity  of  the  time  interval  to  changes  in  laser  aspect  angle 
(a)  is 


a (At) 
da 


Now 
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dx 


(arcsin  u)  = 


du 
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so  that 


( tan  r2 
\ tan  0 


tan  7 2 

tan  a 
da 


(52) 
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The  sensitivity  of  the  laser  aspect  angle  to  roll  rate  is 


& tan  A I A sin2a  . 1/2  . 

5; ' rCI  r 2 ' ~~ur + 11  cosai  • 1 

sin  a + ijtan  71  q 

1 , 

Sensitivity  of  the  laser  aspect  angle  to  time  interval  variation  is 
given  by 


d o _ 


-2ntotan  y 


1 I A sin  a . 1/2 


d(At)  2 . 2 |^7F+T^  cos0t 

sin  a + qtan  y ^ | q 


The  sensitivity  of  laser  aspect  angle  to  angular  spacing  between  the 
retroreflectors  (f3)  is 

do  _ 2 tan  yl  \ A sin2a  1/2 

dp  .22  l/2  + ^ cosa  • (£ 

sm  a + r]tan  7^  q 


Combining  Equations  (54)  and  (56)  through  (60),  one  obtains 


2 / sina  sec2  7^  \2 

(d  - cr)  = ( 5 5 

Vsin  a + qtan  7 j 


~j  ^/2  " “^2  ^cosa+A^2 


• ~ 2 “1/2  \ 2 
sina  sec  72^  • \ 

~ 2 - 2 J 

sm  a + qtan  7 / 


2 2 

(A  + cosa)  (67  ) 


2nAt  tan  7 \ /.  . 2 n /0  \2  „ 

_1_»  / A sm  q , 1/2  \ ,2 

2 2)1  2/2  Tl  cosa)  (5w) 

sin  a + qtan  7^/  \ q / 


(2«w  tan  7. 

”1 I 

sm  a + qtan  7^ 


A sin2a  . 1/2  _\2  ,ca.  v2 

i/J”  + ^ cosal  (SAt) 

4 / 


2 tan  71  \ (A  sin2 a 1/2  \2  2 

. 2 7 “2  \ 172“  11  coscv  (8P)  • 

sm  a + qtan  y^f  \ q j 
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d.  Conclusions 

Although  a complete  system  study  for  the  two  station  con- 
cept has  not  been  completed  the  studies  conducted  to  date  indicate  that 
this  approach  may  be  a cost  effective  range  instrumentation  technique. 
Additional  studies,  similar  to  those  described  for  the  three  station 
system,  are  warranted.  It  should  be  noted  that  if  the  technique  actually 
works  as  described  here  all  of  the  advantages  originally  envisioned 
with  the  use  of  the  three  station  system  would  accrue  to  the  two  station 
system. 
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